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In the Standard Model of particle physics, lepton flavour and lepton number are conserved 
quantities although no fundamental symmetry demands their conservation. I present recent results 
of searches for lepton flavour and lepton number violating hadron decays measured at the B factories 
and LHCb. 

In addition, the LHCb collaboration has recently performed a search for the lepton flavour violating 
decay t~ — > u~ u~ u + . The obtained upper exclusion limit, that has been presented in this talk for 
the first time, is of the same order of magnitude as those observed at the B factories. This is the 
first search for a lepton flavour violating r decay at a hadron collider. 



I. INTRODUCTION 

In the Standard Model both lepton number as well 
as lepton flavour are conserved quantities [U [2] ■ Since 
both can be broken in extensions of the standard model, 
observation of either of them would be a clear sign for 
new physics. 

Results of the search for lepton number violation 
(LNV) and lepton flavour violation (LFV) in decays of 
hadrons are presented. These comprise the B decay 
modes B+ -> h + £ + f-, B+ -> hrt + l' + , and the corre- 
sponding D modes D+ -> h + i + l'~ and D+ -> hrl + i' + . 
The final state meson h may hereby either be a stable 
meson (pions or kaons) or in case of B decays also a D 
meson. Additionally to these modes, the first limit on 
the branching fraction B + — > D°7r + /i~/i~ is presented as 
well as a new result on the search for LFV in the decay 
of t leptons at the LHC. 

Throughout this document charge conjugate decays 
are implied. 

A. Lepton Number Violation 

Numerous models without lepton number conservation 
have been proposed, see [3] for an overview. Similar to 
the fundamental diagram in the neutrinoless double beta 
decay, any neutrinoless hadron decay with two same sign 
leptons in the final state probes the existence of Majorana 
neutrinos. Of the two lowest order diagrams for LNV in 
meson decays, one can go through an on-shell neutrino, 
while the other contains a virtual neutrino (Fig. [T]) . In 
B + decays one of them is Cabbibo favoured depending 
on the final state. Thus for Majorana neutrino masses 
in the accessible mass range (up to 5140 MeV/c 2 ) the 
modes B+ -> T>-£+£+ and B+ -> D*~£ + £+ arc more 
sensitive and also provide a mass measurement. Beyond 
the accessible mass range other final states (n~£ + £ + , 
T>~£ + 1 + ) are more sensitive. 

In the framework of LNV through a fourth neutrino 
iV with a large Majorana mass, an observation of LNV 
not only provides information about the mass 777,4 of the 
fourth neutrino, but also on the WN£ coupling strength 




(a) on-shell Majorana neutrino 




(b) virtual Majorana neutrino 



FIG. 1: Lowest order diagrams for LNV in meson 
decays, involving on-shell or virtual Majorana neutrinos. 
Depending on the individual quark flavours either of 
them can be Cabbibo favoured. Reproduced from [4]. 



I Va I . A compilation of different exclusion limits is shown 
in Fig. [2| For the coupling to the muon the strongest 
constraints come from kaon physics. 

Complementary to the modes with one meson in the 
final state, it has been suggested in [5] to also consider 
B + — > D°7r + /i _ /i _ with the diagram shown in Fig. [3] 
Until 2012 no limit on the branching fraction of this 
decay has been measured. 
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FIG. 2: Constraints on charged lepton couplings V e 4, and 
Vu4 to a fourth heavy Majorana neutrino from 2009 as a 
function of the mass 7714 [3] • 




FIG. 3: Tree level Feynman diagram for the decay 
B" 4 D°7r+/x-/x-. 
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FIG. 4: Feynman diagram for lepton flavour violating 
meson decays in the Standard Model with neutrino 
oscillation. 



B. Lepton Flavour Violation 

In contrast to violation of lepton number, lepton flavour 
violation has been observed experimentally in the neu- 
trino sector. Through loop diagrams, neutrino oscillation 
can also enter the charged sector as illustrated in Fig. [I] 
- the predicted rates however are immeasurable small, 
suppressed by powers of m^/m^r [5]. 

Two examples how to introduce sizeable lepton flavour 
violation are multi Higgs extensions by means of new 
scalar particles (see e. g. [7J) as in the diagram in Fig. 5a 
or by means of heavy neutrinos as introduced in low scale 
seesaw models (e.g. [8]) which couple to electrons and 
muons as shown in Fig. |5b| Other ways to embed LFV 
in the standard model are given e. g. in 0. 

Particularly interesting about LFV in B decays com- 
pared to D decays is that the B mass is high enough 
to produce a r + , [iT pair in the final state. For new 
physics introduced in a Higgs coupling, this final state is 
most sensitive due to the high masses, and thereby Higgs 
couplings of the leptons involved. 

Similarly the decay t~ — > is not entirely for- 

bidden, but neutrino oscillation at loop level alone cannot 
bring the branching fraction to an observable level. As 
presented in [TU] the strongest limits on LFV in lepton 
decays come from the — > e~j mode, the search for 
r~ — > is particularly interesting because some 

new physics models (e.g. Littlest Higgs [H]), as in Fig. [6j 
have enhanced lepton flavour violating couplings to heavy 
leptons (favouring r over \i decays) and do not involve 
photon couplings and therefore enhance the three lep- 
ton final state over the t'y final state [12]. Moreover to 
identify the character of new physics, a search in both 
I -> £'j and £ -t £'£'£' must be performed. 



II. EXPERIMENTAL RESULTS 

The study of rare decays naturally needs large event 
samples, which, for B and D mesons, is available at the B 
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(b) Standard Model extension with a heavy neutrino vh 

FIG. 5: Examples for introduction of lepton flavour 
violating in meson decays 




(a) standard model with neutrino oscillation 




(b) new physics with a charged Higgs 

FIG. 6: Feynman diagrams for t~ — > pT fx~ in 
different models. 



factories and at the LHC. The most stringent constraints 
on LFV and LNV in modes involving electrons come from 
BaBar and Belle, while muonic final states are now best 
constrained by recent LHCb measurements. 



A. Limits on Lepton Number Violation 

Decays of B + , D + , and K + mesons were used to search 
for Majorana neutrinos of different masses. The mass 
difference of the decaying meson and the final state lepton 
is the upper limit on the mass of the on-shell neutrino 
which can be produced. Since the neutrino mass is the 
invariant mass of the final state meson-lepton pair, the 
sum of their rest masses is the lower limit on the accessible 
mass range. 

The strongest limits on the lepton coupling \V e 4\ 2 
and | Vfj,4 1 2 to a fourth neutrino are in the low neu- 
trino mass region between 140 MeV/c 2 and 353 MeV/c 2 
coming from searches for the decays K + — » c + e + 7r~ 
and K + — » /i + /j + 7t~ respectively. Couplings down to 
\Vei\ ;$ 10 -8 are thereby ruled out in the most sensitive 
range. 

The currently most stringent limits on LNV in charm 
decays and thereby higher neutrino masses have been ob- 
tained by the BaBar collaboration (T3] shown in Tab. [Taj 

The extension of the search range to higher masses 
is only possible in B decays, the enormous production 
cross section in hadron collisions makes the LHC the 
optimal place for searches for LNV in B decays. LHCb 
recently provided new results on the on-shell modes B + — > 
7r~u + u + and B + — » D7u + u + , as well as the virtual 



modes B+ -s- D+^V 4 " and B+ ->■ D*"/x+/x+ @]. 

Limits on the branching fraction are hereby set as a 
function of the neutrino mass for the on-shell modes. For 
comparison, assuming a flat distribution of the decay 
products in phase space, the observed branching fraction 
is shown in Tab. [lb] along with the modes which are 
sensitive to virtual Majorana neutrinos and previous 
measurements. 

The first search for B + — > D°7r~ fi + fi + has been per- 
formed by LHCb [3] and showed no excess over the back- 
ground. Since this channel involves an on-shell Majorana 
neutrino, the limit is given as a function of the neutrino 
mass as well. 

LHCb also provides the strongest limits on |V^4| up to 
the B + mass consideringthese results come from B + — > 
7r~^i + ^ + , shown in Fig. [7} 

A natural way to search for a lepton number violating 
decay is to search for same sign leptons from a common 
vertex. Thereby the analysis' implications on an inter- 
mediate on-shell neutrino are only drawn correctly if the 
lifetime of the neutrino is short enough not to degrade 
the reconstruction. To estimate how the observed limits 
are to be understood in models with long lived heavy 
neutrinos, LHCb also provides the relative reconstruction 
efficiency as a function of the neutrino lifetime, shown in 
Fig. |5| 
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TABLE I: Current limits on lepton number violating charm (a) and bottom (b) meson decays. 



(a) charm decays (b) bottom decays 
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FIG. 7: Limit on | V^, 4 1 2 from B~ -> n + n~ 'fi~ measured by LHCb g]. 
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FIG. 8: Relative reconstruction efficiency as a function of the Majorana neutrino lifetime. The branching fraction 
limits from [4] have been computed for the assumption of infinitively short lifetimes (100% relative efficiency). For 
longer lifetimes, the reconstruction efficiency decreases and the observed limit has to be scaled down. The efficiencies 
are given for (a) B+ -> 7r~^" l >+, (b) B+ D7^+/i+, and (c) B+ -> D°7r _ ^ + /i + . 
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B. Limits on Lepton Flavour Violation 

The tightest constraints on lepton flavour violating 
processes in charm decays have been found by the BaBar 
collaboration, listed in Tab. |IIa| For bottom decays, 



Tab. lib shows the recent results, involving r leptons in 



the final state. Details are given in [TS]. These results 
improved the limit on the energy scale at which LFV can 
occur [Tj|| significantly. 

The implication for new physics is that the energy scale 
for LFV effective operators is pushed up from 2.2 TeV 
to II TeV or from 2.6 TeV to 15 TeV for the b ^ d and 
the b — >• s transition respectively [19] . 

The most recent result in this talk is the limit on LFV 
in t~ — » /j," n~ fi + achieved by LHCb [22]. The hadron 
collider environment introduces special experimental chal- 
lenges compared to the B factories. 

a. t tag At the B factories, r are produced in pairs. 
A clean event selection therefore is to look at events with 
four tracks - three from the signal candidate and one 
from a standard model one prong r decay. At the LHC 
the main source for r is the leptonic D~ — > t~v t decay 

eg. 

b. Normalisation The r tag automatically provides 
the number of produced t which enter the analysis. Since 
the main production mode for t at the LHC does not 
provide any further charged particles, the number of r 
entering the analysis is not directly accessible. A normal- 
isation to allowed r decays is not possible since they are 
indistinguishable from more abundant D + decays with 
7T° in the decay chain. 

c. Background Having no production tag, back- 
ground from events without r, such as B and D cascade 
decays, is more severe in the LHCb analysis than for the 
B factories. 

The main advantage of LHCb however is the huge 
production cross section for r from D s decays. Consid- 
ering the charm and bottom production cross sections 
measured by LHCb [23] [24] and the known semileptonic 
branching fractions [35], about 8 x 10 10 r leptons have 
been produced at LHCb in 2011 compared to a total of 
10 9 r pairs at the B factories. 

The analysis strategy of [22] is similar to other rare 
decay searches at LHCb. A loose cut based selection 
is applied to get a processable data sample. All events 
passing this selection are classified in a three dimen- 
sional likelihood space. The discriminating variables are 
the invariant mass of the t~ — > candidate, a 

multivariate classifier A^3body for the three body decay 
properties (geometry, displacement, track quality, isola- 
tion, and kinematics), and a multivariate classifier for 
the particle identification Mpid (combining information 
from muon stations, RICH detectors, and the calorimeter 
signature). The latter classifiers use boosted decision 
trees [26j with adaptive boosting [27] as implemented by 
TMVA HE]. 

The signal efficiency of the multivariate classifiers 
as well as the invariant mass resolution come from 



simulation and are calibrated on a control channel - 
D+ — » 7r + 0(/i + /i~) in the case for the three body classi- 
fier and the invariant mass and B + — > J/iJj(ji + /i~)K + for 
the particle identification. 

The D+ — > 4>tt + calibration channel also serves as a nor- 
malisation, since the branching fractions B(D+ — > 4>tt + ) 
and B(D+ — > t + v t ) are known - yielding the number 
of t which have been produced in D s decays. To deter- 
mine the fraction of r from D s decays, /(D s ), the bottom 
and charm cross sections measured by LHCb [531 HI] i 
as well as the branching fractions of charm and bottom 
hadrons to r are used. Hereby most of the systematic un- 
certainties (e. g. luminosity measurement, reconstruction 
efficiencies) cancel, i.e. /(D s ) is more accurately known 
than the inclusive r production cross section. Contri- 
butions from gauge bosons or Drell-Yan processes have 
been evaluated to be negligible. 

Using the above normalisation as well as the efficiencies 
for selection, reconstruction, and trigger the branching 
fraction can be written as follows: 



A 4 M M J = 



6( D+ -> <KmV> + ) 



x /(D 8 ) 



N T - 



l/i-t-/i-j7T-t 



where e r , 



and select the normalisation decay and e s i g 



is the total efficiency to trigger, reconstruct 

is the total 

efficiency for the signal channel 

The dimuon decay of the <f> is chosen to provide similar 
trigger and particle identification properties compared to 
the signal being sought for. The non resonant contribu- 
tion from D+ — > /i + /i~7r + decays was found to be below 
2%. 

The three dimensional likelihood space is subdivided 
into 150 bins (five for .M3body and A^pi p, an d six for the 
invariant mass) as shown in Fig. [9] and 10a The signal 
efficiency for each bin is evaluated from the calibration 
channels and the background in each bin is estimated 
from the sidebands in the invariant mass. 

The background consists mainly of two components. 
Firstly combinatorial background which is modelled by 
an exponential and secondly by D+ — > ?7(/i + /i _ 7)/i + i^ 
decays. This physical background is not discriminated in 
the current analysis by either A^3body or Mpin as it has 
the same behaviour as the signal in all input quantities. 
Rejecting this decay will be subject of future improve- 
ments. It is modelled by an exponential multiplied by a 
second order polynomial for which all shape parameters 
have been fixed on simulated events. The normalisation 
is left free in the final fit within one standard devia- 
tion from the expected yield which is determined using 
the normalisation channel and the branching fractions 
B(D| -> W +^),B(p -> M + /^7), and B(D+ -+ <f>n+). 
For illustration the invariant mass distribution and the 
combined fit for the combination of the two highest 
■M3body and two highest VWpid bins is shown in Fig. [TOb 
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TABLE II: Limits on lepton flavour violating hadron decays at 90 % confidence level. All listed limits from the BaBar 

collaboration. 
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(b) bottom decays 
limit 
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(a) Distribution for simulated background and the simulated 
signal as a function of the 3 body decay classifier. 
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(b) Distribution for simulated background and the simulated 
signal as a function of the PID classifier. 



FIG. 9: Distribution of signal events in the two multivariate likelihoods for signal (blue / solid) and background (red 

/ dashed). 



TABLE III: Limits on the branching fraction for 
r~ — > fi~ fi~ [i + obtained by different experiments. 



collaboration 



limit 



Belle < 2.1 x 10" 8 @90% CL [3T] 

BaBar < 3.3 x 10" 8 @90 % CL [32] 

LHCb < 6.3 x 10~ 8 @90% CL [22] 



For the final limit, all bins are combined using the CLs 
method [531 130] ■ The observed limits at 90 % confidence 
level is 6.3 x 10~ 8 , in agreement with the exp ected limit 
for the absence of a signal (8.2 x 10~ 8 ). Tab. Ill shows 



III. CONCLUSION 

Hadron decays measured at the B factories and at the 
LHC provide an excellent and abundant probe to search 
for LNV and LFV. So far no signal has been observed 
and only lower limits for the branching fractions are 
given. The B factories have achieved high sensitivity and 
ruled out branching fractions to the level of 10 -5 . The 
huge production cross section of B mesons at the LHC 
furthermore enabled LHCb to improve the limits on LNV 
in B decays further. Finally, the first LFV r decay search 
performed at a hadron collider has been performed. 



the comparison to the limits from BaBar and Belle. 
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(a) simulated signal (b) observed events 



FIG. 10: Invariant mass distribution for (a) simulated signal candidates and (b) observed events in the two highest 
■^3bod y and A^pid bins. The background fit (D+ — » r]fi + u in green / dotted; combinatorial in red / dashed; 
combined in blue / solid) is shown in the range which is used for the fit. 
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